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1. Introduction 

In SU (N c ) vector gauge theories, perturbation theory allows for a non-trivial infrared fixed 
point (IRFP)[jlJ, ||]. When an IRFP exists in the very strong coupling region, investigation by lattice 
gauge theories is needed for going beyond perturbation theory. There have been a lot of recent 
lattice studies; see Refs. for reviews. 

We are motivated by Ref. []|] in which a systematic lattice analysis over a wide range of the 
gauge coupling constant was performed. According to those authors, Nf = 3 in SU (2) and Nf = 7 
in SU (3) belong to the conformal window. Their result was based on the prediction of the pion 
mass behaviour in the strong coupling limit, given in Ref. ^ as follows: 

. i (l-16fc 2 )(l-4>c 2 ) 
C ° Sh ^ = 1+ 8^(1-6^) ' (U) 



where K = 2m * +8 and mq is the bare mass parameter of Wilson fermions. Eq. ( |1 . 1| ) implies that the 



pion becomes massless at K = 1/4, independent of N c and Nf. Therefore the authors of Ref. [Q] 
investigated at j8 = ^ = in particular, which led to their prediction for the critical flavour number 
of the conformal window mentioned above. 

However, with Wilson fermions the phase structure in the strong coupling limit is very compli- 
cated, as pointed out in Ref. due to the existence of a parity-flavour broken phase (Aoki phase). 
In the Aoki phase K > k c = 1/4 at j3 = 0, and on the phase boundary the charged pion is massless 
while the neutral pion is not 1 . 

Therefore, in Ref. [0] and this proceedings we check the particle spectroscopy in the strong 
coupling limit to determine whether Ref. [||] is correct, and we explore the phase structure of the 
Wilson fermion system with many flavours to determine whether the Aoki phase exists in the many 
flavour case. 



2. Simulation details 

We use Nf flavours of mass-degenerate Wilson fermions [f7]]. Dynamical configurations are 
generated by the standard Hybrid Monte-Carlo (HMC) algorithm with Nmd^t = 1. Lattice sizes are 
6 2 x 12 , 8 x 16 , 12 x 24 2 and 12 3 x 24 though not all of these sizes were used for every choice 
of Nf', this range of sizes allows us to look for finite size effects (for some Nf). In this exploration, 
we compute observables from 50 ~ 100 configurations with 4 ~ 5 intervals between trajectories 
after thermalizing. We monitor the following observables: m 2 K , m p , the plaquette value((D)), the 
axial Ward-Takahashi identity quark mass (m^ WI = ^^^^^w^ )» the Polyakov loop, Creutz 
ratios, the chiral condensate (or the propagator norm), the lowest eigenvalues, and (S(t)S(O)) as a 
function of t. Due to page constraints in these proceedings we focus on a few observables most 
relevant to a discussion of the phase structure. 

3. N c = 3 case: Comparison with the data of Ref. [|5p 

Iwasaki et al.[|5|] concluded that the Nf = 6 case is in the confinement region because they 
found Ncg > 10000 in the Molecular-Dynamics (MD) evolution for the thermalizing process (they 

'The order parameter of the Aoki phase is (ij/y 5 T3y/). 
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did not attain the thermalization) and that the plaquette value was decreased toward zero, in contrast 
with the case of Nf > 7. Thus they interpreted these observations to be signals for a massless pion, 
in spite of no direct calculations of the pion mass and the quark mass. However, as plotted in the 
left panel of Fig. [T], we found results for the pion and quark masses as functions of 1 jx at j8 = 
that differ from Ref. [§]: we did attain thermalization with N CG < 10000 in MD step at K = 0.25 
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Figure 1: 51/ (3) gauge theory with Nf flavours at j3 = 0. Left panel: m| and m^ WI as functions of 1/jc. 
The analytic prediction is shown as a solid blue curve reaching m\ = at 1 / K = 4. Right panel: Plaquette 
history in HMC and R-algorithm in various lattice setups. 

for Nf = 6, and our pion mass is not zero as shown in the left panel of Fig. [I]. In addition, our pion 
and quark masses depend on Nf such that only the quenched case (Nf = 0) obeys Eq. ( |l.li) leading 
to the conventional Aoki phase. 

To understand the origin of this discrepancy between our data and theirs, we monitor the 
history of the plaquette value on various lattices with various setups. In addition to our own HMC 
simulations with periodic boundary conditions reported in Ref. |J7j|, we now monitor the plaquette 
value by using MILC code[||]. The MILC code uses the R-algorithm with anti-periodic boundary 
conditions, which are the same algorithm and boundaries as were used by Iwasaki et al. in Ref. 1^]. 
We generate the configurations with a cold start for Nf = 8, then switch to Nf = 7 from the last 
configuration of Nf = 8 and switch to Nf = 6 from that of Nf = 7, and we plot the plaquette value 
history in the right panel of Fig. [j] For larger A^,, the simulation is not affected by the different 
algorithms and boundary conditions. However, we found that the plaquette value at TV, = 4, where 
Iwasaki et al. obtained their conclusion, shows a large deviation. 

Our observations at ./V, = 4 are consistent with the results of Ref. [gj], i.e. the small value of the 
plaquette and the pion mass. However, it seems clear that results with this small time extent is not 
representative of larger lattices; there seems to be some kind of the thermal phase transition. We 
conclude that our studies with larger N t , which showed a massive pion, are more appropriate for 
zero-temperature discussions, while conclusions from Ref. ^ about the conformal window cannot 
be supported. 



4. N c = 2 case: Exploration of the phase structure 

For SU (2) gauge theory, Ref. [Q] predicted that Nf = 2 is in the confinement region and that 
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Nf = 3 is in the conformal window for the same reason as discussed for SU (3): Ncg > 10000. 
However, as we have already pointed out in the previous section, finding Ncg > 10000 during the 
thermalizing process is not reliable evidence for discovering a conformal window. Although we 
are unable to obtain direct results for Nf = 2 due to the very small value of the pion mass around 
K = 0.25, we did obtain results for Nf = 4, 6, 8, • • • and we plot those in the upper left panel of 
Fig. H| for various Nf flavours and in the upper right panel of Fig. ^| for Nf = 0,6 and 12 especially 
to check finite size effects. In the upper panels of Fig. @, there are a confinement phase and an 
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Figure 2: SI/ (2) gauge theory with Nf flavours at j8 = on a 6 2 x 12 2 lattice. Upper left panel: m 2 n and 
mf l as functions of 1/jc. Upper right panel: Volume independence for Nf = 0,6 and 12. The solid blue 

curve represents the prediction of the equation. Bottom panel: jX — yj 'Xo(H^), where Ao(#^) is the lowest 
eigenvalue of H™, as a function of 1/fC. 



unknown massive pion phase, and these two phases are separated by a transition at K = iq 2 . 

Does an Aoki phase exist in the strong coupling region? The SU (2) case is qualitatively the 
same as the SU (3) case: the pion and the quark mass reveal the 2-state signal, the quenched case 



(Nf = 0) agrees with Eq. ( |1.1| ) but the dynamical case (Nf > 0) depends on Nf, and further there is 
no K c beyond which m\ = 0. Namely our result is 



0.25 = Ke(N f = 0) > . . . > Kc(N f = 6) > . . . > K c (N f = 12) > . . . 



(4.1) 



if there exists in Nf > case. Therefore it seems there is no standard Aoki phase due to non 
existence of the massless pion over the whole region of K. 

Further we monitor the lowest eigenvalue in the bottom panel of Fig. ||, }X - 



M H w)> where 



H\v is the hermitian Wilson-Dirac operator, in order to consider the phase structure of many Wilson 

2 This is not a precise value. Hysteresis is observed within a range of K for Nf > 6 in SU (2), and is a label for 
the location of that transition. 
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fermions. The lowest eigenvalue for Nf > is not zero in both the confinement and the unknown- 
massive-pion phase, while the eigenvalue in Nf = tends to go to zero as K" — > 0.25. According to 
the modified Banks-Casher relation^], the order parameter of the Aoki phase is (\(f(x)Y5^3Y(x)) = 
2npH w {h = 0) where pji w (X) is the eigenvalue distribution function of H w . The data from small 
lattices suggest that the order parameter of the Aoki phase is zero, and this is consistent with the 
absence of a standard Aoki phase in the case with many Wilson fermions. However, it would be 
necessary to investigate on the larger lattices and to calculate the spectrum of low-lying eigenvalues 
before making a strong conclusion, and this is left to future studies. 



Q— O In confinement 
A -A In unknown phase 
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Figure 3: m\ as a function of not only positive but negative K for Nf = 6 at j3 = on 6 2 x 12 2 . 

The outline of the phase structure in our research is shown in Fig. ||. The existence of a massive 
pion phase (the red triangles in the figure) is apparent, but without that information one might try 
to obtain fQ. by extrapolating to m 2 K = 0. In that way one would find that there is the gap of K c 
(fc+ ^ K~) for m 2 = between the positive and the negative region. According to the traditional 
interpretation, this is the Aoki phase containing a massless charged pion. However, the massive 
pion phase (red triangles) indicates that, at least with the many Wilson fermions (Nf/N c S> 1), 
another phase has appeared instead of a standard Aoki phase. 

Therefore, in Ref. [J7|] and this proceedings, we don't have a clear signal for the existence of 
an Aoki phase and the Sharpe-Singleton first order scenario [[H]], but instead we have obtained 
evidence for a different phase at strong coupling. 

5. Summary and Discussion 

We explored SU (2) and SU (3) gauge theories with many Wilson fermions and reported result 
in Ref. [|7|] and additional results in these proceedings. In contrast to Ref. [g], we find that (i) The 
pion mass in the confinement region depends on Nf. (ii) As defined by extrapolation to m 2 = 0, 
Kc depends on Nf and shifts from 0.25. (iii) No massless pion exists for any value of K. Thus 
the extrapolated K c is not a very useful quantity, and extrapolation to m K = is not valid when 
searching for an Aoki phase, (iv) There is the 2-state signal and hysteresis around K = JQ, which 
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depends on Nf. (v) The Nf = 6 case in SU (3) gauge theory with Wilson fermions is not in the 
confinement region because the pion mass at K c = 0.25 is not zero. Simulations at N t = 4 are not 
representative of larger N t values, so we do not validate the claim[]|] that the critical flavour of 
the conformal window is 7 (Nf ^ 7). (vi) There is a standard Aoki phase in the quenched case 
(Nf = 0), but the situation is markedly different in the large flavour case (Nf/N c 3> 1). 

Now we make an effort to identify the unknown massive pion phase. The Polyakov loop is 
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Figure 4: Left panel: Polyakov loop normalized by the quenched case, (|L|)/\y — (\L\)^. = q, as a function 
of 1 / K. Here we take the absolute value of the Polyakov loop. Right panel: Creutz ratio, J ( 1 , 1 ) , with 1 
APE smearings as a function of 1/fC. 

plotted in the left panel of Fig. [|. Along the small extent (N x = 6 on 6 2 x 12 2 ), the Polyakov loop 
shows finite temperature behaviour with a 2-state signal and hysteresis. However, the Polyakov 
loop along the large extent (N t = 12 on 6 2 x 12 2 ) doesn't show any transitions and is consistent 
with the quenched case (Nf = 0). This means there is no finite temperature transition on large 
lattices, so the unknown phase is not the result of a confinement-deconfinement phase transition. 

The right panel of Fig. [| is a Creutz ratio, #(1,1). Atj6=0 the static quark is confined (modulo 
string breaking), so the Creutz ratio at short distances is near unity in lattice units. We compute 
#(1,1) for K > JQ, and find that it becomes small in contrast with its value in the confinement 
phase. This relative difference might imply a transition from confinement to weak coupling. 

Although we monitored various quantities to identify the phase, unfortunately at present we 
cannot determine the true nature of this unknown massive-pion phase. Still, this phase is the confir- 



mation of a (first order) bulk phase transition at the strong coupling, as pointed out in Refs. [ 11 , |12[ ] 
by observing the high- and the low-plaquette phase and the chiral condensate with the first order 
transition. Our result is the first observation of that longer list of quantities (the pion mass, the 
quark mass and so on) and the fact that the standard Aoki phase seems to be invisible due to this 
bulk phase transition at strong coupling. Furthermore, it is found that this bulk phase transition 
even appears in the strong coupling limit (j8 = 0) with Wilson fermions. This implies that the (first 
order) bulk phase transition is caused by fermionic dynamics only, not by the gauge dynamics, 
namely this phase is independent of the type of lattice gauge actions. 

Although it's difficult to determine which phase is realized, it will be possible to determine 
whether the standard Aoki phase has appeared or not, by observing the order parameter of the 
Aoki phase, (xj/ysTiy) . Refs. Jl3| , 14] carried out a calculation of the order parameter for Nf = 2 
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in SU(3) (Nf/N c < 1) by using the twisted mass Wilson fermions. Let us remember again the 
bottom panel of Fig. |[ When this figure is turned upside-down, a result similar to the data of the 
order parameter might be obtained in SU (2), due to the very rough relation (1//75T3I//') ~ jj on finite 
lattices. If this is true, SU (2) twisted mass fermions might show a bulk phase transition in Nf = 2, 
instead of Aoki phase. Therefore we are calculating the order parameter in Nf = 2 and Nf > 6 in 
SU(2) gauge theories (Nf/N c > 1) with twisted mass Wilson fermions. It should be possible to 
clarify the phase structure of a system with (many) Wilson fermions in the near future. 
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